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We show some recent HERWIG results related to the top quark mass reconstruction at 
the LHC and discuss possible improvements for studies of the bottom quark fragmenta- 
tion function in the top decay. 
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Qs ' The large amount of tt pairs which are expected at the future experiments at the 

f^ , LHC"'^ will allow precise measurements of the top-quark properties. For the sake of 

^^ ' an improved measurement of the top mass rrit , a new approach has been recently 

suggested^, which consists of using the invariant-mass distributions of J/iJj + £ 

i-^ I pairs, where the J/ip is produced by the decay of a 6-flavoured hadron and i is the 

lepton coming from the W produced in the top decay t — > bW. Even though the 

Mh. branching ratio for the process B -^ J/tp is small, it has been shown that it provides 

(— I I a clean signal at the LHC and about 10"^ well-reconstructed final states are expected 

in one year of high luminosity. The mji^i distributions can then be compared with 

a template of shapes parametrized by the top mass and mt can be fitted. The 

PS ' expected experimental uncertainty, which is dominated by the statistical error, has 

^ . been estimated to be Amf 2± 1 GeV. This is a method which crucially relies on the 

Monte Carlo description of top production and decay and on the model which is 

used to simulate the fe-quark hadronization. 

Standard Monte Carlo algorithms'^'^ simulate multiparton radiation in the soft 
or coUinear approximation and leave regions of the phase space completely empty 
('dead zones'). In particular, HERWIG simulates the top decay in the top rest 
frame. We do not have any soft-gluon radiation from the top quark, while the h 
quark is allowed to radiate gluons in the 'forward' hemisphere Q < 6g < 7r/2, with 
9g being the emission angle of the soft gluon with respect to the h direction^. The 
following parton cascade satisfies the prescription of angular ordering^. The total 
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energy loss due to gluon radiation is roughly correct, as the collinear singularity 
0g ~ is correctly taken into account; nevertheless problems are expected in the 
angular distributions, since the W 'backward' hemisphere tt/2 < 6g < tt is a dead 
zone for the shower. 

The new version HERWIG 6.1^ allows gluon radiation in the backward hemi- 
sphere via the implementation of matrix- element corrections*: the dead zone is 
populated according to the 0{as) matrix element ('hard correction') and the par- 
ton cascade in the already-populated region is corrected by the use of the exact 
amplitude any time an emission is the 'hardest so far' ('soft correction')^. 

We wish to use HERWIG 6.1 to study invariant-mass distributions which are 
relevant to extract rrit at the LHC^"'^^. For simplicity, we study the spectra iriBe, 
which must be convoluted with the B — > J/ip decay spectra to obtain the rrij/^i 
distributions. In Fig. 1 we show the iriBe distributions at the LHC according to 
HERWIG before (6.0) and after (6.1) matrix- element corrections to top decays, for 
mt — 175 GeV; we observe a shift towards lower values of mse after the inclusion 
of hard and large-angle gluon radiation. In Table 1 we show results for the average 
values {mBt) for different values of mt and observe that the shift induced by matrix- 
element corrections is of about 800-900 MeV. If we parametrize the relation between 
{niBe) and mt according to a straight line, we find, by means of the least-square 
method: 

6.1 : {mBi) = 0.568 mt- 6.004 GeV , e(GeV) = 0.057 ; (1) 

6.0 : (mBi) = 0.559 mt - 3.499 GeV , e(GcV) = 0.052 , (2) 

e being the mean square deviation in the fit. 

Inverting the above relations, we find that the values of mt extracted using the 
two versions HERWIG 6.0 and 6.1 differ by about 1.5 GeV, a value which is larger 
than the expected uncertainty. This implies that, for an accurate reconstruction 
of the top mass, the corrections to the top decay must be applied. If we set a cut 
m.B£ > 50 GeV to reduce the effect of possible backgrounds on the low-mass tails, 
the impact of matrix-element corrections to top decays is reduced to 1 GeV, but it 
is still competitive with respect to the expected uncertainty^^. 

Table 1. Results for the invariant mass nigi at the LHC for different values of mt- 



mt 


(ms.)"'^ 


a(6.1) 


(ms.)"-" 


cr(6.0) 


(ms;)"'" - {msi)"-' 


171 GeV 


91.13 GeV 


26.57 GeV 


92.02 GeV 


26.24 GeV 


(0.891 ± 0.038) GeV 


173 GeV 


92.42 GeV 


26.90 GeV 


93.26 GeV 


26.59 GeV 


(0.844 ± 0.038) GeV 


175 GeV 


93.54 GeV 


27.29 GeV 


94.38 GeV 


27.02 GeV 


(0.843 ± 0.039) GeV 


177 GeV 


94.61 GeV 


27.66 GeV 


95.46 GeV 


27.33 GeV 


(0.855 ± 0.039) GeV 


179 GeV 


95.72 GeV 


28.04 GeV 


96.51 GeV 


27.67 GeV 


(0.792 ± 0.040) GeV 



The next step of the presented analysis is a more detailed study of the b frag- 
mentation in the top decay ^^. Another possible approach to obtain the mBi distri- 
butions is to perform an exact next-to-leading order calculation of the decay rate 
t -^ bW{g) and convolute the result with the b fragmentation function, taken from 
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Fig. 1. mgi distributions at the LHC according to HERWIG 6.1 (solid histogram) and 6.0 (dotted). 

some LEP or SLD data"'^^. It will be very interesting to compare the HERWIG 
results to the ones of this calculation. In order for such a comparison to be trust- 
worthy, we have to tune the parameters of the HERWIG cluster model to fit the 
e+e^ data. We have a recent tuning by the OPAL collaboration^'', after which the 
agreement of HERWIG to the LEP data is considerably improved. However, even 
after such a tuning, problems have been found in the comparison to the SLD data^^, 
which was not taken into account as an input for the OPAL tuning. The tuning of 
the HERWIG hadronization model using both LEP and SLD data is in progress^^. 
In summary, we reviewed some HERWIG results which are relevant to the top 
mass determination at the LHC using final states with leptons and J /^ and found 
that the effect of matrix-element corrections to top decays is Awt ~ 1.5 GeV. We 
also reported on further studies for a deeper understanding of the h fragmentation 
in the top decay, using data from the e^e~ machines. 
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